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Newsletter No. 10 



Category III Trigonometry Modules 



February 1, 1971 



We are enclosing four trigonometry modules. The module "CIRCULAR 
FUNCTIONS" was written by Mr. Robert Gillespie of Allderdice. Its place 
in mathematics is described in the overall plan he and his group prepared 
last summer. The "12th grade mathematics" topics given in Newsletter No. 

7 came out of that plan. 

4 

As a result of mathematical questions raised by "CIRCULAR FUNC- 
TIONS", two other modules were prepared called "TRIGONOMETRY FUNC- 
TIONS AND TCHEBYCHEV POLYNOMIALS" and "INVERSE CIRCULAR 
FUNCTIONS". The module "VORTAC" is for the more ambitious' student. 
Students who want to try a two aircraft simixlation of navigating by VORTAC 
should type: 

>RUN 166FW /WGAME/ 



Call For Student Programs 

We would like to circulate student programs of all levels throu^ the 
newsletter. This is not a request for "best programs", but rather a mech> 
anism through which teachers can give recognition to experienced as well 
as neophyte student-users of confuting. Please include a listing, sample 
run, and appropriate documentation, including teacher’s name, subject, 
student’s name, age, and grade. 






ERIC 

M/lUffllilRlffTlM 



0- 

CO 

o 

o 

(J 



Exchange Programs 

We are looking into the possibility of exchanging newsletters and pro- 
grams with other high schoo’. computing efforts. The two newsletters we are 
aware of right now come from the Southern Minnesota School District Project 
at Mankato, Minneapolis, and from Brooklyn Polytech in New York (The 
Huntington Project). This later group is specializing in simulation (type II) 
programs. 



^Supported in part by IKF grant GJ 1077 
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CIItCULAR FUNCTIONS 



Consider an angle A with its initial side placed along the 
X-axis, with its vertex at the origin F:(0,0), and with a terminal 
side that passes through the point F:(X,Y). Let us call the dis- 
tance from (0,0) to (X,Y) R, so that = x2 + Y^* Then the six 

circular functions of the angle A are defined as: 



Sine of the angle = y/r 
Cosine of the angle - x/r 
Tangent of the angle = y/x 



Con tangent of the emgle s x/y 
Secant of the angle = r/x 
Cosecant of the angle = r/y 



Sample problem: 



Find the circulcu: functions of an angle (A) whose terminal 
side passes through the point F: (5^12). 



Solution: 


X = 


5 and y 


therefore 


r ss 


SQRT (5*5 


thus 


sin 

cos 


A = 12/13 
A - 5/13 




tan 


A = 12/5 




cot 


A = 5/12 




sec 

CSC 


A = 13/5 
A = 13/12 



* 12 

+ 12*12) » SQRT(169) 

- .9231 

- .3846 

- 2.400 

* .4167 

= *2,600 
= 1.0833 




Using the Computer 



On the next page we will show the flow chart for a program 
that asks the user to supply values for X and Y. The program 
then calculates SIN A, COS A, and TAN A. 

Notice that care is taken to avoid dividing by zero. (For 
which circular functions can this happen???) 

More advanced programming problems are suggested on pages 
5 and 6. 



*N0TE: When using the computer you should type 5 as 5. or 5.0 

and 12 as 12. or 12.0- etc. 



3 



1 



CIRCULAR FUNCTIONS 



Flowchairt for a Program to Calculate Sin, Cos, and Tan 
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Here is an NBS program based on the flow chart on page 3: 

50 PR, "THIS PROGRAM CALCULATES THE SIN, COS, AND TAN FUNCTIONS'* 

60 PR, "OP AN ANGLE A DEFINED BY THE RATIOS Y/R, X/R, AND Y/X ** 

70 PR, "RESPECTIVELY" 

95 PR, 

100 PR, "TYPE IN THE COORDINATES OF YOUR POINT IN X,Y ORDER, USING" 
110 PR, "DECIMAL FORM" 

120 INPUT X,Y 
200 I^SQRT (X*X+Y*Y) 

300 S»Y/R 
400 C=X/R 

450 IF X=0 GOTO 900 
500 T=Y/X 

600 PR, "SIN(A)=**:Sj"COS(A)=";C;"TAN(A)=":T 
700 PR, "DO YOU WISH TO ENTER ANOTHER POINT"? 

710 INPUT R$ 

720 Y$="YES, SURE, OK, O.K,, CERTAINLY, YUP,OF COURSE, RIGHT, YEAH, Y, " 

730 IF IEQIV(R$,Y$,1) GOTO 100 
800 STOP 

900 PR, "SIN(A)=":S?"COS(A)=":C?"TAH UNDEFINED" 

910 GOTO 700 
>RUN 

THIS PROGRAM CALCULATES THE SIN, COS, AND TAN FUNCTIONS 
OF AN ANGLE A DEFINED BY THE RATIOS Y/R, X/R, AND Y/X 
RESPECTIVELY, 

TYPE IN THE COORDINATES OF YOUR POINT IN X,Y ORDER, USING 
DECIMAL FORM, 



|>?5,0, 12,0 



SIN (A) - , 9 2 31 COS ( A) = , 3 84 6 TAN ( A) -2 , 4 

DO YOU WISH TO ENTER ANOTHER POINT? SURE 
TYPE IN ETC, 



‘* NOTE : You must type in two decimal numbers- If you don’t know a 

number (say the square root of 3) , you can use @NBS, 

EXAMPLE : Suppose you want to find the circular functions for 
P:(/3T , /2 ), 

>RUN 



WRONG!! 



TYPE IN THE COORDINATES OF YOUR POINT 
>?SQRT(3), SORT (2) 



>RUN 

THIS PROGRAM ETC. 



CORRECT 



TYPE IN THE COORDINATES OF YOUR POINT IN X,Y ORDER, 
?@NBS 

>PR, SORT (3) 

1,732 

>PR, SORT (2) 

1,414 



>EXIT 

^?1,732, 1.414 
* ETC 



0 




1 
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Computer' Problems 



1. Use tlie program on Page 4 for P:(l,l), P: (0,1), P:(1,0), 

P:(-3,3), P:(-2,3), P: (-3,"3) , P:(3,-3) £IMPT: Type 1 as 1.0 etc.] 

2. Modify the program given on Page 4 so that all six circular 
functions are calculated. 

a. Use your computer program to find the circular functions of 
the following: 

i. An angle A whose terminal side passes through 
P: (2*SQRT(2) , SQRT(7)> NOTE: You must use §NBS 

to find 2*SQRT(2) etc. 

ii. An angle A whose terminal side passes through 
Ps(-2.386, 7.590) 

iii. An angle A whose terminal side passes through P: (0,2.7) 

iv. An angle A whose terminal side passes through Ps (2.7,0) 

b. Run your program again, supplying data on several points in 
all four quadrants. Can you make a general rule about the 
signs of the circular functions in each quadrant? 

c. Investigate angles whose terminal sides are coincident with 
the coordinate axes; that is, circular functions of 0, 90, 

180, 360". Which functions are undefined? 



Advanced Problems : 

3. Ifrite a program that automatically generates the six circular 
functions 

FOR X ~ -1.0, .9, .8, .7, «6, «5, «4, «3, .2, <1^ 0. 

*.l, *.2, *.3, “.4, *.5, *.6, “.7, *.8, *.9, *1.0 

WITH R FIJffiD AT R = 1 (UNIT CIRCLE) . 

Here is part of a program to do this with some output: 

100 AS="2(SD.DDB) 2<SD.DDDB) 3<SDD.DDDB) SDD.DDDB/'* 

no BSs'»2<SD.DDB> <>(SD«DDDB) * UNDEF • SDD.DDDB • UNDEF ' SDD.DDDB/** 
120 CSs'*2(SD.DDB) 2(SD«DDDB) SDD.DDDB ' UNDEF * SDD.DDDB* UNDEF */** 
130 Rsl 

140 PRINT ** X Y SIN C0S TAN C0T SEC CSC** 

150 F0R X=1 BY -.1 WHILE X>-1.1 
160 IF ABS(X)>sl G0T0 280 
170 Y=ABS<SQRT<R+R-X*X)) 

180 IF ABSCXX.OOS G0T0 250 
190 IF ABSCYX.OOS G0T0 280 

200 PRINT IN F0RM ASiXl Yl Y/RiX/Rl Y/XlX/YiR/XiR/Y 
210 NEXT X 
230 END 

o 
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3. (Continued) 


























>RUN 
















X 


Y 


SIN 


C0S 


TAN 


COT 


SEC 


CSC 


+ 1*00 


+0*00 


+0*000 


+1 *000 


+00*000 


UNDEF 


+01*000 


UNDEF 


+ 0*90 


+0*44 


+0*436 


+0*900 


+00*484 


+02*065 


+01*111 


+02*294 


+ 0*80 


+ 0.60 


+0*600 


+0*800 


+00*750 


+01*333 


+01*250 


+01*667 


+ 0*70 


+0*71 


+0*714 


+0*700 


+01*020 


+00*980 


+01*429 


+01*400 


+ 0*60 


+ 0*80 


+0*800 


+0*600 


+01.333 


+00*750 


+01*667 


+01 *250 


+ 0*50 


+ 0*87' 


+0*866 


+0*500 


+01.732 


+00*577 


+02*000 


+01 .155 


+ 0*40 


+0*92 


+0*917 


+0*400 


+02.29J 


+00*436 


+02*500 


+01*09 1 


+ 0*30 


+ 0*95 


+0*954 


+0*300 


+03*180 


+00*314 


+03*333 


+01 .*^48 


+ 0*20 


+ 0*98 


+0*980 


+0*200 


+04*899 


+00*204 


+05*000 


+01*021 


+ 0* 10 


+0*99 


+0*995 


+0* 100 


+09*950 


+00* 101 


+10*000 


+01*005 


+ 0*00 


+ 1*00 


+1*000 


+0*000 


UNDEF 


+00*000 


UNDEF 


+01*000 


-0* 10 


+ 0*99 


+0*995 


-0* 100 


-09*950 


-00*101 


-10*000 


+01*005 



ETC* 



4 . Mem Against Machine : 

Try to make the computer blow its mind by calculating 
TAN (A), with R = 1, and X approaching zero* 

10 LET X=1 

20 IF X<«0000000000001 G0T0 60 
ONE 30 PRINT ••X=”lXlTAB<35>fTAN<X)=”jS0RTn-X*X>/X 

APPROACH 40 LET X=X/2 
50 C0T0 20 
60 END 



5* The picture on the cover shows that the paint dripping from 
th^ paint disi^nsers on the tliree wheels (rotated by monkey 
power) traces out "graphs'* of the circular functions on the 
canvas (moved by mad artist power) . 



Write a 




(a) 


OR 


(b) 


OR 


(c) 



program that: 



Graphs the SIN or COS function* 

Graphs the SIN and COS function on the same axis, using 
two different values of R* 

Graphs SIN and COS functions with different values of R 
and different starting points, possibly making "computer 
generated" art as your goal* 



Suggestions : 

1. Ask your teacher for the module called "Graphing by Computer* 

Could' you make your graph print symbols like "2", "3", etc* 
to make a "number" painting for later water coloring, where 
1 = Red, 2 = Blue, etc* 



2 * 



INVERSE CIRCULAR 





But the principal value of ARCSIN(0.8) is 
0,9273 RADIANS WHICH CAN BE CAL CULATED BY 

• ARCSIN(X) TT/2 - Ul - X" ) * 
(i.570795207- .214512362 * X 
+ .087876311 * X*- ,044958884 * X’ 

+ .019349939 * X*- .004337769'* X') 
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FINDING THE INVERSE CIRCULAR FUNCTIONS 



In the module on trigonometric functions we saw that 
polynomial approximations cure used by computers to calculate 
the circulcir (or trigonometric) functions when the indepen- 
dent variable is the angle. The simple NBS program: 



>10 LET A=.9273 

>20 PRINT ''A=":A;"SIN(A)=*':SIN(A) 

>RUN 

will result in the output: 

A^.9273 SIN(A)=. 800002869 

In order that the computer successfully execute this program 
it must use one of the Hastings approximations for finding 
SIN(.9273). (See Problem 3 in the module "TRIG. FUNCTIONS 
AND TCHEBYCHEV APPROXIMATIONS".) 

THE INVERSE PROBLEM 

Now suppose we ask the computer to 

>PR. ARCSIN(.8) 

we will obtain: 

0.927295218 

We have asked the computer to print: 

"the angle which has. a sine of 0.8" 
or to be exact: 

"Print the angle A (in radians) such 
that SIN (A) —0.8, and A lies in the 
interval -ir/2<A£ir/2 

The cover of this module shows why we have to agree on an 
interval for A. The inverse function ARCS IN is not s inale 
valued, and so we must eliminate confusion by using the 
(unspoken) assumption that the pyincipal value of A will be 
calculated. For the SINE function, ^e principal value is 
found in either quadrant I or quadrant IV. 



HASTINGS SHEET 35 



The formula for approximating ARCS IN <X) shown on the 
front cover is Hastings' Sheet 37, a rather accurate for- 
mula* It assumes that X is limited to the range 

0<X<1 

Here is an NBS program to use a simpler but less ac- 
curate formula (Hastings' Sheet 35) to calculate the ARCSIK 
function* The program also prints the more accurate NBS 
ARCSIN function for compcirison* The Sheet 35 formula is: 

ARCSIN (X)^/2- tm [) * (1* 5707288- * 2X21144*X+* 0742610*X*X 
-*0187293*X*X*X) 

Notice that polynomial is evaluated using nested multiplica 
tion in line 560 in order to reduce the nui^er of multipli- 
cations* Why is this inportant? 



PROGRAM' .0 COMPARE NBS LIBRARY ARCSIN FUNCTION 
WITH HASTINGS 35 FUNCTION 



argument NBS LIBRARY HASTINGS 35 



0 

0 * 1 
0*2 
0*3 
0*4 
0*5 
0*6 
0*7 
0*8 
0*9 
I 

>LISTNH 

500 PR*'*PROGRAM TO COMPARE NBS LIBRARY ARCSIN FUNCTION '* 

501 PR*'*WITH HASTINGS *35 FUNCTION" PR* PR* 

510 LET A0*l *5707288 LET Als-*2I2II44 

520 LET A2s*07426l0 LET A3*-*0I87293 
530 P=3* 14159/2 

5 40 PR*"ARGUMENr”:TAB<20) f'NBS L I BRARY": TAB < 40 ): "HASTINGS 35” 

545 PR*” ” 

550 FOR X«0 TO I STEP 0*1 
5 60 D=<<A3tX+A2)4(X+AI)+X + A0 
570 S=P-S0RTC1-X)*D 

580 PR*X:TAB(20)$ARCSIN(X):TAB(40) :S 
590 NEXT X 
600 END 



0*100167421 
0*201357921 
0*304692654 
0*41 1516846 
0*523598776 
0*643501109 
0*775397497 
0*927295218 
1 * 1 19769515 
1 *570791655 



6*620000204E-05 
0* 100107027 
0*201313803 
0*304700575 
0.41 1559864 
0*523643617 
0*643521213 
0.775385477 
N 0*927265853 
1*119751877 
1.570789605 



1 



4 




Pr Ob lean 1 

Improve the above program so that X may lie in the 
interval -l£X£l, Test your program by running it with line 
550 changed to: 

550 FOR X=^l TO 1 STEP 0,1 
[HINT : RECALL THAT SIN (-X) =-SIN (X) ] 

Problem 2 

Find out what NBS does when you type 
>PR, ARCSXN(2) 

Can you make your program do the same for the Hastings 35 
approximation? Or do you have to? 



Problem 3 (OPTIONAL) 

Change your program to use the Hastings Sheet 36 for- 
mula: 



ARCS IN (X)*tt/ 2- (>fT^)*(l, 57078786-0. 21412453*X+0, 08466649*X*X 
-0 , 03575663*X*X*X+0 , 00864884*X*X*X*X) 



Use 'fT=3, 14159265 



Problem 4 (OPTIONAL) 

Run the program for Problem 3 again, using Hastings’ 37 (on 
cover) and the Double Precision feature of NBS, (See NBS Primer, 
Section 7, "EXTENDED DATA TAPES’*, Use ir=3 , 141592653589793238 ) 



Problem 5 



In practice, the library functions of NBS are used for 
calculation. Practice using the ARCSIN(X), the ARCCOS (X) , 
and the ATAN(X) functions*in NBS to find the angles A, B, C, 
D, E, F in the figure below. (Use direct mode if you wish.) 






% 






Trigonometric Functions & 
Tchebychev Approximations 

(OR — CAN A HUMBLE COMPUTER THAT ONLY KNOWS HOW TO 
DO ARITHMETIC CALCULATE TRANS C^DENTAL FUNCTIONS?) 



Problem 1. Below you see three representations of the "sine" func- 
tionT Which is the ri^t one? 

□ 1. (a) 

□ 2. (b) . 

□ 3. (c) 

□ 4. All Of the above. 

□ 5. None of the above. 



(a) The SINE 

function : 



(b) The SINE 

function : 



(c) The SINE 

function: 



Y 


Y/R 


X 


±y(R"-X")/R 


A 


SIN(A) 


0.0 


0.000 


1.0 


0.000 


0.0 


0.000 


C.l 


0.100 


0.9 


0.436 


0.1 


0.099 


0.2 


0.200 


0.8 


0.600 


- . 0.2 


0.198 


0.3 


0.300 


0.7 


0.714 


0.3 


0.295 


0.4 


0.400 


0.6 


0.800 


0.4 


0.389 



etc . etc . etc . 






Note: We are using the 

letters X, Y, A, and R 
to represent the same 
quantities as in the mo- 
dule "Circular Functions". 
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• * 

* TCHEBYCHEV 

Answer to Problem 1; All of the representations are correct. 

Here are three NBS programs which you can try to convince your- 
self that you can get three different- looking ”gre^hs” of the sine 
function. The shape of these graphs depends on which variable (X, 
y, or A) you decide to increment in equal steps. 

After looking over these programs, you may get the feeling that 
something has been *put over* on you. Question: Has a little mathe- 

matical "hanky-panky" crept into these demonstration programs???? 
Which one??? (For the answer, read onlM) 



Demonstration Program (a) 



2 






+2 
+ 2 
+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 2 
2 + 

+ 



2 



2 



200 
210 
220 
2 30 
2 AO 
2 50 
2 60 
270 
2 80 
290 

295 

296 
2 99 



F0R Y*0 STEP *1 UNTIL Y>*95 
PRINT TABOO) TA8(30+30:*=Y)s "2" 

NEXT Y 

F0R Y=1 STEP UNTIL Y<*05 

PRINT TAB(30) J TAB( 30+30+Y) :*'2" 

NEXT Y 

F0R Y=0 STEP -*1 UNTIL Y<-*95 
PRINT TABCSO+SO^Y) :*'2":TAB(30)?'+*' 
NEXT Y 

F0R Y=-l STEP *1 UNTIL Y>*05 
PRINT TAB(30+30*Y) :"2":TAB(30) :•'+*' 
NEXT Y 
END 



+ 

2 + 






* 
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Demonstration Program (b) 




+ 1 




+ 


1 


+ 


1 


+ 


1 


+ 


1 


+ . 


1 




1 


+ 


1 


+ 


1 


+ 


1 


+ 


1 


+ 


1 


+ 


1 


+ 


1 


A 


1 


+ 


1 


+ 


1 


+ 


1 


+ 


1 


+ 


1 


1 + 


_ 


1 + 




1 + 




1 + 




1 + 




1 + 




1 + 


100 F0R X=1 STEP -*l UNTIL X<--95 


1 + 


no LET SsSORT< l-X*X) 


+ 


IRO PRINT TABOO) :”+”:TAB<3b+30*S):”l” 


+ 


130 NEXT X 


> 


140 F0R X=-l STEP •! UNTIL X>-95 


+ 


150 LET S^-SORTU-XtX) 


+ 


160 PRINT TAB<30+30tS)r*M”jTABO0):'"<” 


1 + 


170 NEXT X 


1 + 
■ ^ 


180 END 



1 + 
I + 

I + 

1 + 



14 
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Demonstration Program (c) 




300 F0R A*0 STEP 6.283/40 UNTIL A>6.283 
310 LET S=SIN(A) 

320 IF S>0 THEN PR. TAB( 30> TAB( 30+30*5) : ”3” ELSE PR. 
TA8< 30+30*5) : ”3”: TAB( 30) * ”+” 

330 NEXT A 
340 END 



+ 

+ 

+ 3 

+ 3 

3 + 

3 + 

3 + 

3 + 

3 + 

3 + 

3 + 

3 + 

3 + 

3 + 

3 + 

3 + 

3 + 

3 + 

3 + 

3 + 

3 + 



3 



3 



3 




3 



3 + 

3 + 

3+ 



3 



3 



3 



3 



3 



3 



If you haven't spotted the haaiky-panky by now, check line 310 
above. Assuming computers can only do arithmetic (which is true), 
an esqpression like 30 + 30 * S in line 320 makes sense, but it 
isn't clear that arithmetic is being used in line 310* Where is 
the eurithmetic being done? Answer: SECRETLYl The computer actually 
goes into a "library" sub^routine at line 310, which tells it what 
arithmetic operation must be done on the angle A (esqpressed in 
radians) in order to calculate SIN (A) . 



n CO n 
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The Secret Unveiled 



The library functions in digital computers used to evaluate 
complicated functions almost always employ expertly designed poly- 
nomials as approximations to these functions. These approximations 
can be made quite accurate, using techniques originated by the 
Russian mathematician Tchebychev (also spelled Chebyshev) , and put 
into useful form by Hoistings.^ 

Here is one of the approximations given by Hastings: 

(Hastings Sheet 14) 

Function: sin 

Range: -1 < X < 1 




Approximation : 




« C^X + C3X* 




where = 1.5706268, C 3 = -.6432292, and = .0727102. 



Before using this approximation, it will be useful to make two 
changes. First, we will change the range of the independent vari^le 
by the substitution: 

y - (ir/2) * X so that ^dien: -1 < X < 1 

we wij.1 have: -ir/2 < Y < ir/2 



Note: Both Y and X are angles , measured in radians. 

. Don't confuse these variables with the Y and X 
coordinates used on the cover. The variable 
Y in the rest of this discussion is identical 
to the variable A shown on the cover. 



After making this substitution, the Hastings Sheet 14 approxi- 
mation becomes: 

sin(Y) = Kj^y + 

vrtiere Kj^ - (2/ir) * Cj^i K3 - * C3; Kg = (2®/ir®) * Cg. 



^A book your library ought to have is Hastings, Cecil, Approximati ons 
for Digital Computers , Princeton University Press, 1955. 
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The second change we will make is to rearrange the right side 
of the preceding formula into "nested multiplication" form: 

sin(y) (((Kg * y * y + K 3 ) * y) * y + k^) * y 

The advantage of this form is that it only takes five multiplica- 
tions, whereas the original polynomial takes nine. 

Here is an NBS program which uses this form to calculate 
sin(y), and compares it with the NBS library routine (which uses a 
higher degree approximation — see problem 3) . 



PROGRAM T0 COMPARE NBS LIBRARY SINE FUNCTION WITH 
HASTINGS I4« 



RADIANS 


NBS LIBRARY 


HASTINGS 14 


-1 


-0.841470985 


-0.841534239 


-0.9 


-0.78332691 


-0.783406887 


-0.8 


-0.717356091 


-0.717433028 


-0.7 


-0.644217687 


-0.644277689 


-0.6 


-0.564642473 


-0.564678888 


-0.5 


-0.479425539 


-0.47943851 


-0.4 


-0^389418342 


-0.389413182 


-0.3 


-0.295520207 


-0.295505151 


-0.2 


-0. 198669331 


-0. 19865316 


-0. 1 


-0.099833417 


-0.099823323 


0 


0 


0 


0.1 


0.099833417 


0.099823323 


0.2 


0.198669331 


0.19865316 . 


0.3 


0.295520207 


0.295505151 


0.'4 


0.389418342 


0.389413182 


0.5 


0.479425539 


0.479438509 


0.6 


0.564642473 


0.564678888 


0.7 


0.644217687 


0.644277689 


0.8 


0.717356091 


0.717433028 


0.9 


0.78332691 


0.783406887 


1 


0.841470985 


0.841534239 


>LISTNH 






290 PR.' 


"PROGRAM TO C0Mt»ARE NBS LIBRARY 


SINE FUNCTION WITH 


291 PR.' 


"HASTINGS 14.** PR. PR. 




300 LET 


P=3. 14159265 




310 LET 


P2sP*P LET P3*P*P2 




320 LET 


Cls(1.5706268*2)/P 




330 LET 


C3=< -.6432292*8 )/P3 




340 LET 


C5»<. 07271 02*32)/<P2+P3) 




350 PR. 

DP 


**RA01ANS*'ITAB(20)|""NBS LIBRARY' 


"|TAB( 40) t *"HASTIN6S 



360 FOR Y*-l TO I STEP 0.1 

365 IF ABS(Y)<*000l LET Y^O 

370 SI=<<<C5*Y*Y+C3)*Y)*Y+CI)*Y 

380 PR. •' "lYJTAB<20)tSIN<Y)lTAB<40)lSI 

390 NEXT Y 

400 END 



1 




ERIC 

MliffllilBIffTlTEB 



0 

y 
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Problem 2. Change the FOR loop in the above program to cover the 
range for Y = 1 to 3 radians# and see what happens ^ Also try 
the range for Y = -1 to -3. If sometliing "blows up,” can you 
modify your program so that any value of Y can be ^ed in the 
Hastings' approximation? (Hint: recall sin(Y) = sin(Tr-Y) 

= sin(2TT+Y) = sin(3TT-Y) = etc.) 

- Problem 3. Here is the Hastings Sheet 15 approximation ^diich 

ou^t to come a lot closer to the values produced by the NBS 
library function. Modify your program to try it out# again 
compcuring it with the NBS library function. 

sin((TT/2) * X) * X + C 3 * + Cg * 

where 

C 3 - 1.570794852 C 3 - -0.645920978 

C 3 = 0.079487663 = -0.004362476 

-1 < X < 1 

Probl^ 4. (Optional) If you find the idea of super-precision of 
interest # ask your teacher to see the Com-Shcu:e XTRAN Library 
Manual# and esqplore the use of double precision curithmetic 
and double precision library functions. One of the best hand- 
books to use as a standard for “correct" values is "The Hand- 
book of Mathematical Functions# Applied Math Series No. 55" 
available from the Superintendent of Documents# U.S. Government 
Printing Office# Washington# D. c. 20402 ($6.50). * A paperback 
version published by Dover Co. may be available at book stores 
for about $5.00. 

Probl^ 5. (Easy) Use the NBS library functions TAN and SQRT to 
find the lengths of the guy-wires X and Y shown below. (Use 
of direct mode NBS as a desk calculator is a good way to 
solve this problem.) • / / 




/ 





Phantom VORTAC 

^ This module will guide you in preparing the master program 
for an on-board flight computer. The computer takes informa- 
tion from a VHF omni range system (VORTAC ) , and calculates the 
magnetic course and distance to a given airport for the pilot. 
The output of the program is such that the pilot can fly to- 
ward a "phantom** VORTAC located at any airport he selects. 

^ A description of this newly developed navigational system 
and the mathematics on . which it is based are given on page 2. 
^ Pages 3 snd 4 outline; alternate methods for handling the 
computation/ and suggest ways in which previous programs you 
have written might be incorporated as sub-routines. 

^ A "real time" simulation of a flight using this system is 
suggested as an advanced level program. 
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Pilots flying over the United States (and most other countries 
of the world) rely on radio facilities called VORTACs for navi- 
gational information. The basic information the pilot receives in 
the cockpit is his position relative to the VORTAC, given in polar 
coordinates . 

The pilot in the illustration below would describe his position 
(obtained from his radio instruments) as being, "on the 100** radial 
of the Allegheny County (A6C) VOR, 50 miles out". 

It is easy for this pilot to 
note that he can get to A6C by 
turning right, and flying a course 
of 280® (Why 280®?)* If he is 
going 200 miles per hour, it is 210 
also easy for him to estimate 
that he will arrive at A6C in 
15 minutes (Is this exactly true? 

— go back over the module on vec- 
tor addition if you are not sure.) 

The catch to all of this is 
that the location of the VOR is 
usually different from the location of the destination airport. 

This difficulty can be handled by flying tfro legs, the first from 
the present position to the VOR^ and the second from. the VOR to the 
airport* This is obviously an inefficient route. 




V 




A new navigational system 
uses an on-board, special-purpose 
computer to tell the pilot what 
course and distance to fly in or- 
der to go directly from his present 
position to the airport. Let's 
first examine three ways of analyzing 
the mathematics involved in such a 
computation . 



Very high frequency Opii Rwge and TACA N* where TACAN is an older 
mTlita^ system of distance measuring equipment. Most civilian 
pilots call these facilities VOR-DHE stations. 

**The angular direction of the intended fli?^t path, measured clock- 



wise from N. 
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Before writing a program for such an on-bOcu:d computer, it 
would be useful for a programmer to solve by hand typical problems 

At this point, let's look at three such 



arising in this situation, 
solutions. 



D=210® 



VORTAC 

E=135® 




Problei^ 1 . The aircraft is 
150 miles out on the 135® radial 
hnd the destination airport is 
100 miles out on the 210® radial. 
That is, in the diagram. Si = 150, 

E = 135®, S3 = 100, D = 210®. 

Let A2 be the angle determined 
by Si and S3. Hence, A2 = 

210® - i35® = 75®. By the Law of 
Cosines, 



S2 = /lOO^ + 150^ - 2 



100 



150 • cos(75®) = 157.3 miles. 

How by the Law of Sines, sin a 3 » (100 * .966)/157.3 » .614 and 
cos A3 = /l -.614^= .789. From this we can get a 3 = cu:ctan(.614/.789) 
= 38®. (Why wasn't A3 computed directly from sin A3?) Therefore 
X ss 135® + 180® - 38® = 277®. Output to the pilot is 157.3 miles, 
277®— is this answer reeisonable? 

Problem 2 . The airplane is 250 miles out on the 45® radial 
and the destination airport is 100 miles out on the 280® radial. 

Ih teacms of the diagram. Si = 250, E = 45®, S3 = 100, D « 280®. 

Converting the positions of the aircraft and destination to 
rectangulcu: coordinates with the VORTAC as origin we have 

sin 45® = 176.8 miles, 
cos 45® = 176.8 miles, 
sin 280®« -98.48 miles, 
cos 280®« 17.36 miles. 



X - 250 
P 

y„ = 250 
P 

= 100 



100 



A = arctan 



176.8 - 17.36 
176.8 - (-98.48) 
arct^Ji (.5792) = 30.1® 




VORTAC 



This is another way of saying that computes do not remove the re- 
sponsibility of analyzing a problem before solving it — quite the 
contrary; they demand, more thought than ever. This may be, in fact, 
one of the most important contributions computing systems can make 
to leaaming. 
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Therefore X = 270® - 30 . 1® = 239 .9® . By the Pythagorean Theorem 
S2 = /(176.8^ 17.36)2 + (176.8 - (-98. 48))* ^ 318.1. Hence the 
output to the pilot is: 318.1 miles, 239.9®. 

Problem 3. The aircraft is 150 miles 
out on the 120® radial and the destination 
airport is 400 miles out the 150® ra- 
dial. **In the diagram, B - 120®, Si = 150 
miles, D - 150®, S3 - 400 miles. Using 
the same reeisoning as in Problem 1 yields: 

A2 = 30® 

S2 := /400* + 150*- 2 • 150 • 400 • cos 30® 

= 280.3 miles 

sin A3 = 400 . (.5236/280.3) « .7135 
cos A3 = /I - .7135* = .7007 
A3 = arctan (.7135/. 7007) = 45.52® 

X « 120® + 180® - 45.52® = 254.48® 

Output to the pilot is: 280.3 miles, 

254.48® — does this answer look reason- 
able? What went wrong? 

Assignment. Write a program for an 

on-bOcu:d navigational computer which will 

accept cis input the radial and distance from a VORTAC of both an 

aircraft and a destination, and which will compute a course and 

distance to the destination. The trigonometric subroutines in the 

computer require arguments in radians, but pilots think in terms 

of degrees, so it will be necessary for you to convert degrees to 

radians and back (see the module on converting to radians) . Other 

possibly Tiseful topics are inverse trigonometric functions , the 

* 

Laws of Sines and Cosines , and transformation of polcu: to rectan- 
gular coordinates . . 



* 

See modules you have worked previously. 

** 

By now you should have noticed that the word radial is used to 
designate the angular position of a line segment that starts at the 
VORTAC. 



N 

1 
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Optional Section. During most flights, the aircraft moves 
through the air and the air moves as well. Write an addition to 
your program which will: 

1. Accept as additional input, specified by the 'pilot', 

(a) the aircraft ''s speed, (b) the speed of the wind, (c) the direc- 
tion of the wind, (d) an elapsed time, t, and (e) a heading * for the 
aircraft. Any heading should be acceptable: the 'pilot' should be 
able to fly wherever he likes. 

2 . Compute the new position of the aircraft on the basis 
of the above information. 

3. Repeat steps 1 and 2 cis often as desired or necessary to 
reach the destination airport. 

YOU can look at the input/output sequence of one program that 
acconplishes this and has been written and stored. 

In order to gain access to this program, pleeise logon (direc- 
tions are given in an introductory module) using the User Identifi- 
cation Code given to you by your teacher: 




f: 



In order to interact with this stored program in the statements: 



Prom this point on the operation of the program should be self- 
explanatory. 



* 



Heading is defined as the angular direction of the longitudinal 
axis of the aircraft with respect to North. In the picture below, 
the pilot is flying a course of 90* , but his heading is 65*. Why? 
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GROSS FLOW CHART FOR SUGGESTED ASSIGNMENTS 



Aircraft 

Position 

Airport 

Position 




